In silico design of metal-dielectric nanocomposites for solar energy applications 
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Recently, a homogenization procedure has been proposed, based on the tight lower bounds of the 
Bergman-Milton formulation, and successfully applied to dilute ternary nanocomposites to predict 
optical data without using any fitting parameters [Garcia et al. Phys. Rev. B, 75, 045439 (2007)]. 
The procedure has been extended and applied to predict the absorption coefficient of a quaternary 
nanocomposite consisting of Cu, Ag, and Au nanospheres embedded in a Si02 host matrix. Sig- 
nificant enhancement of the absorption coefficient is observed over the spectral range 350-800 nm. 
The magnitude of this enhancement can be controlled by varying the nanosphere diameter and the 
individual metal volume fraction with respect to the host matrix. We have determined the optimal 
composition resulting in enhanced broadband (350nm-800nm) absorption of the solar spectrum us- 
ing a simulated anneaUng algorithm. Fabricating such composite materials with a desired optical 
absorption has potential appUcations in solai^ energy hai^vesting. 

I. INTRODUCTION 

Judicious selection of the materials used for a given photonic application is of paramount importance 
in tailoring the desired optical properties of a composite material [1]. Light waves, when directed to the 
interface between a metal and a dielectric are capable of a resonant interaction with the mobile electrons 
present at the surface of the metal, giving rise to surface plasmons, i.e. density waves of electrons that 
propagate along the interface 121]. Similarly, such nanocomposites offer much potential in designing ma- 
terials with enhanced light absorption due to the local electromagnetic field enhancement from the surface 
plasmon resonances that occur at the metal-dielectric interface jsl]. The wavelength-dependent absorption 
of such metal-dielectric nanocomposites can be tuned by varying parameters such as the shape, size and 
volume fraction of the nanoparticles as well as the dielectric constants of the host and metal |l4i bL lal- 
The increasing emphasis on renewable energy sources has generated interest in finding cost-effective ways 
of enhancing the efficiency of solar energy harvesting systems. One promising route is the use of metal 
nanoparticles or nanocomposite coatings as an absorbing layer on silicon devices, thereby giving rise to 



photocurrent enhancement that correlates well with the plasmonic response of the applied coating 116, 



tzlli 



While consideration of the effects of electromagnetic scattering due to the coating morphology on the de- 
vice performance needs to be addressed [8], optimization of the absorption of incident solar energy though 
design of an absorption coating is a critical step in the development of improved solar absorbing mate- 
rials. Recently, we have developed and validated against experimental data a hierarchical homogenization 
procedure capable of accurately predicting the absorption characteristics of multi-metal/dielectric nanocom- 
posites 1 9]. Typically, optical analysis and experimentation in the literature are limited to nanocomposites 



with one or two differing species qLZLIsL 



10, 



nil . Here, we show that by employing the above homogeniza- 
tion procedure each metal added in the dielectric host introduces an added degree of freedom from which 
the optical response of a multi-metal nanocomposite material may be tailored. In particular, we demon- 
strate how a Cu:Ag:Au-Si02 nanocomposite might be used as an absorption coating optimized to the solar 
spectrum via the manipulation of parameters such as volume fraction and particle size. The ability to de- 
sign such materials from first principles can guide future experimental efforts aimed at contro lling pa rticle 
size and spacing in techniques such as laser-induced dewetting Ivl. \l3 , 
laser-assisted chemical vapor decomposition 1 15, 1^ . 
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II. HOMOGENIZATION PROCEDURE AND NANOCOMPOSITE SPECTRAL RESPONSE 



The Bergman-Milton formula lllTl 11811 for the effective dielectric constant for a binary metal-dielectric 
composite is given by: 
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where / corresponds to the component's volume fraction, e is the component's dielectric constant with 
subscripts a and h referencing the metal and host matrix respectively. 7 is a geometrical factor taking 
into account the shape of the metal particles. For the purposes of this paper we are interested in spherical 
particles and accordingly choose 7 = |(1 — /q). Equation [T] denotes the effective dielectric function for 
the two component case of a metal particle embedded in a host dielectric matrix. Note it is customary to 
let //i = 1 — ^ /„ where corresponds to the metal components. In Fig. \T\ the effective absorption 
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-Jm{eeff ) where A is the wavelength of incident electromagnetic energy, of a 



coefficient, a = „ . , 

single metal species of either Ag, Au, or Cu embedded in a dielectric Si02 matrix is presented. The effective 
absorption coefficients of Ag and Au in Si02 have the desirable property of exhibiting strong absorption 
near their plasmonic resonances; although this response is localized at 414 nm and 529 nm for Ag and 
Au nanoparticles respectively. Cu embedded in SiOi exhibits a rather broad spectral response spanning 
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from 200 nm to 645 nm. However, the maximum obtainable magnitude of the effective absorption of 
the Cu-Si02 system is clearly weaker than that of Ag-Si02 and Au-Si02. Quantitatively, over the range of 
wavelengths spanning from 225 to 1200nm, the Cu, Ag, and Au in SiOi systems in Fig. [T]have an integrated 
absorption / = J22^l^!^ (^dX of 0.229, 0.338, and 0.203 (unitless) respectively. Thus, Ag exhibits the 
most enhancement over the spectral range, but this enhancement is localized and may not be the optimum 
response for a given application. 

For our proposed design, we limit the analysis to a dilute system (i.e. total volume fraction of metal 
versus dielectric < 10%) because it is desirable from a practical standpoint of cost and fabrication efficiency. 
A second constraint we have placed is to set the minimum nanosphere diameter to be no less than 10 nm. 
Above this size regime, quantum confinement effects may be neglected II19I1 . Moreover, nanoparticles of 
size > 10 nm can be assembled on Si02 substrates using robust nanomanufacturing techniques such as laser- 



induced dewetting |[l2[[l3 . 14], ion implantation I lo[ 
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llh . and laser-assisted chemical vapor decomposition 



1611 . In Fig. |2j the homogenization procedure is presented for a dilute 3-metal mixture in a host 
dielectric matrix. The top level represents the effective dielectric function for the target nanocomposite of 
design interest. At each level, the total volume fraction is constrained so fh + fn = 1- The schematic 
for homogenization of a 2-metal mixture and alloys in a host dielectric matrix has been presented elsewhere 
iQ]. Here, the 2-metal mixing rule is extended to treat the quaternary system by viewing each component as 
having an effective permittivity (^^ff, ^e^f ^e//- During the mixing process, the average electric field 
within the composite is held constant while the final effective permittivity is calculated using equal volumes 
at each level of mixing. As described in Ref. |^: the effective permittivity of an N-component mixture 
can be determined by mixing N-1 binary mixtures, each comprising of a host and a distinct metal, the host 
being common to the N-1 pairs. 

The nanoscale effect of change in the electron relaxation time due to scattering of electrons on the 
metal-dielectric interface can be quantified via a modified Drude model |20] expressed by Eq. |2l 
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where Tg// is the effective electron relaxation time, Ti^^ik is the bulk electron relaxation time, dn corresponds 
to the particle's diameter and vp is the Fermi velocity of the electrons. The size effect may then be accounted 
for by modifying the imaginary portion of the effective dielectric coefficient [91] as follows 
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where ujp is the plasmon frequency and oj is the angular frequency of the incident electromagnetic energy. 
By exploiting this nanosize effect, one is able to broaden the plasmonic peaks at the expense of their mag- 
nitude as demonstrated in Fig. [3] 

Ag and Au embedded in Si02 exhibit very intense absorptions near their plasmonic resonance. On the 
other hand, the spectral response of Cu is broadband, extending into the range of green light. The net effect 
of mixing these particular materials allows one to achieve a wide degree of freedom in choosing the spectral 
response of the composite. Figj4] shows the nonlinear response of the composite system as a function of 
variations in the volume fraction. Here the volume fraction of each metal is permuted between 6%, 3%, and 
1% while the size of the particles is fixed at 15 nm. Changing volume fraction results in a nonlinear change 
in the wavelength depended absorption. Further, decreasing the diameter of a given metal broadens the 
spectral behavior of the metal's plasmon resonance at the expense of absorption magnitude. Fig. [3] shows 
the trade-off between spectral broadening and plasmon resonance magnitude. To demonstrate the wide 
degree of control of this particular quaternary system. Fig. [5]presents examples of manipulating magnitude 
over the spectral range as a designer might do for a given application. Our mixing approach predicts the 
intrinsic spectral behavior of the nanocomposite with certain volume fraction ratios and particle diameters. 
Extrinsic parameters, such as the actual composite volume and area exposed to solar irradiation, determine 
the magnitude of solar energy absorbed. 



III. OPTIMIZATION 



Towards designing the initial framework of an solar-harvesting device, we have used a simulated an- 



nealing algorithm 1121 



22, 



23h to determine the optimum composite configuration. We define the "en- 



ergy" to be minimized as E{{, d) = ^ abs \{si — ai(f, d))| with composite parameters f =/i, f2--fn and 



d =di,d2--dn. Si and Oj represent the normalized solar spectrum data ll24ll and effective nanocomposite 
absorption at a given wavelength i respectively. Since the majority of the solar energy is concentrated in 
the visible band (see Fig. |6ll, we only define this "energy" for i spanning from 350 to 800 nm. This choice 
of E allows for the absorption profile that most closely resembles the shape of the solar spectrum to be 
determined. The algorithm begins by initializing the system to some initial state E^{i, d). A neighboring 
state E^ is called by using the condition 



^"(r,d") = s^(f + Af,d + Ad) 



(4) 



where Af = fmax'^ with x being a vector of dimension n = 3 with each component taking a random 
value between {—1,1}. fmax represents the maximum magnitude the volume fraction of a particular species 
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may step. Similarly, Ad = dmax ■ x with dmax representing the maximum magnitude the diameter of the 
particles corresponding to a particular species may step. Here, values of fmax = 0.1% and dmax = 0.1 
nm provided sufficient finesse in moving through the search space. Next, the neighboring state is compared 
with the best state encountered thus far E'^ (where E'^ was originally initialized to the same state as E^). If 
j^n ^ j^b jj^gj^ j^b equal to E"-. Then, the algorithm must decide if this neighbor state will become the 
preferred state for the system. To accomplish this, a Boltzmann-type probability P analogous to classical 
statistical physics is calculated- namely 

^ = exp( -(^;-^^h (5) 

where T is the annealing parameter analogous to temperature and A: is a constant used to refine the 
annealing schedule. Note that if P > 1 then P is simply reassigned to 1. P is then compared to a 
random number x between {0, 1}. If x < P then is accepted as the new system state and E^{{, d) = 
P"(f", d"). Notice from Eq. [5]that if E"- < E^, then the neighboring state is always accepted as the new 
system state. The process of calling neighbor states and deciding whether or not to accept them as system 
state is repeated over C cycles. This represents a random walk of C-steps through the {f , d} parameter 
space. For optimization calculations performed here, C was chosen to be 500 as for values greater than 
this no appreciable difference in the results was observed. At each iteration of the algorithm, only values in 
the dilute regime and outside the realm of quantum effects were accepted (i.e. X] /n ^ 0.1 and dn > 10 
nm were enforced at each step). In addition, the maximum particle size was constrained to dn < 30 nm as 
mentioned above and only physically realizable values of volume fractions were permitted (/„ > 0). If the 
algorithm suggested a neighbor state outside the allowable domain, a large value was assigned to the energy 
(E"- ~ 10^). This, by virtue of Eq. |5l gives a very low probability for such states to be accepted. Finally, T 
is reduced and the entire process repeated again. T is reduced according to a prescribed schedule until it is 
nearly equal to zero after iterations. Notice that as T is slowly reduced the system begins to accept lower 
and lower energy configurations until the it is forced to into a (global) minimum. For calculations conducted 
here, the annealing schedule prescribed was T/v = (1— /i)r/v-i where /i = 0.99 using N = 1000 iterations. 
The initial value of T was chosen so that the probability of the algorithm proceeding from a lower state to 
a higher state and vice versa was approximately the same. This ensured that the search space was relatively 
large during the initial stages of the anneal. A value of T = 100 with A; = 0.1 was used for simulations 
here. Typically, f = < 3.0%, 3.0%, 3.0% > and d = < IS.Onm, IS.Onm, IS.Onm >was taken as an 
initial guess. When numerical parameters were determined such that different simulation runs yielded the 
same answer (within a tolerance of AE ~ 10^^), different initial states where chosen and simulations 
ran again. This was done to verify that the algorithm had, in fact, determined the system configuration 



5 



yielding the global minimum. The optimized configuration is presented in Fig. [6] with material parameters 
1.80%Cu:0.35%Ag:6.4%Au with respective particle diameters of 10.1 nm, 29.6 nm and 10.0 nm. 

IV. CONCLUSION 

In summary, we have extended our homogenization procedure to quaternary nanocomposites. This 
procedure yields an analytic expression that acts as a constraint for tuning the spectral characteristics of the 
nanocomposite system. Via manipulation of metal volume fraction and nanoparticle diameter, one is able 
to tune the spectral response over the degree of freedom allowed by the nanocomposites' constituents. We 
note the selection of Ag, Au, and Cu in SiOi allows for a controllable broadband filter with a wide degree 
of freedom over a spectral range of 350-800nm. Furthermore, we have shown how a nanocomposite of this 
type maybe optimized to offer a promising potential application as an absorption coating on solar devices. 
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Figure 1: The wavelength-dependent response of single species of 30 nm diameter nanoparticles embedded in Si02. 
Ag and Au exhibit strong optical absorption at 414 nm and 529 nm respectively. Cu in SiOi yields a fairly broad 
resonance spanning from 200 nm 645 nm. 




Figure 2: Schematic of nanocomposite mixing rule to find effective dielectric constant for the case of dilute quaternary 
metal-dielectric system. 
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Figure 3: Dependence of composite absorption on particle size calculated from homogenization procedure. The 
volume fractions of the metal species are held fixed at 3%. The size of the nanoparticle diameters are permuted 
between 10 nm, 15 nm, and 20 nm. Decreasing the diameter of a given metal broadens the spectral behavior of the 
metal's plasmon resonance at the expense of absorption magnitude. 
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Figure 4: Dependence of absorption on volume fraction calculated from homogenization procedure. The nanoparticles 
for all metal species are held fixed at 15 nm. The volume fraction of each metal species is permuted between 6%, 3%, 
and 1 %. The ensuing nonlinear response of the effective medium's absorption coefficient is shown. 
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Figure 5: Manipulation of broad-band absorption for a dilute Cu:Ag:Au-Si02 mixture. The only constraint is that 
the nanoparticle diameter lie between 10 nm and 30 nm. The spectral response can be tailored to a given application 
fairly precisely over the spectral range of 350-800nm. 




Figure 6: Optimized absorption of the solar spectrum for a Cu:Ag: Au-Si02 nanocomposite. The optimal composition 
of the system was found to be 1.80%Cu:0.35%Ag:6.40%Au with respective particle diameters of 10.1 nm, 29.6 nm 
and 10.0 nm. Spectral data curve fit is plotted in arbitrary units for visualization. 
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